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The distribution of the low-lying QCD Dirac spectrum is analyzed by means of 
partial quenched chiral perturbation theory. We identify an energy scale below 
which the valence quark mass dependence of the QCD partition function is given 
by chiral Random Matrix Theory. This is the the domain where the low-energy 
QCD partition function is dominated by the fluctuations of the mass term of the 
Goldstone bosons. In mesoscopic physics this domain is known the ergodic domain 
and the corresponding energy scale is the Thouless energy. 



1 Introduction 

It has become clear, both from numerous lattice QCD simulations 0i and 
hadronic phenomenology, that QCD at low energies is characterized by chiral 
symmetry breaking and confinement. Because of spontaneous breaking of chi- 
ral symmetry, the low-lying excitations are given by the associated Goldstone 
modes. Because of confinement, these are the only low-lying excitations. In 
this lecture, we wish to point out another important feature of the low energy 
QCD partition function, namely disorder. Ultimately we would like to answer 
the question whether confinement can occur if the classical motion of quarks 
in generic Yang-Mills fields is not chaotic? 

fne way to understand confinement is by analogy with impurity scattcr- 
The argument starts from the semiclassical expression for the quark 
propagator in an external gauge field given by a sum over classical trajecto- 
ries weighted by a phase. Its average over all gauge fields vanishes because 
of this phase. On the other hand, the pion propagator given by the absolute 
value of the quark propagator survives averaging over all gauge field configu- 
rations. We will analyze the quark propagator using a spectral representation 
of the Dirac operator. From the study of disordered systems we know that the 
eigenvalues of localized states with and exponentially falling wave function are 
uncorrelated, whereas eigenvalues of extended states with a power-like fall-off 
of the wavefunction are correlated according to Random Matrix theory. Spec- 
tral correlations therefore constitute a sensitive measure for the tail of the wave 
functions. One of the questions we are interested in is whether the correlations 
of the Dirac eigenvalues support this picture of confinement. 

The spectrum of the Dirac operator and the breaking chiral symmetry are 
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not unrelated. According to the Banks-Casher formula S the order parameter 
of the chiral phase transition is given by 

^ = ^. (1) 

where p(A) is the average spectral density p{X) = {J2k^i^ ~ ^k)), and V is 
the volume of space time. The average (• • •) is over the QCD action. With 
the smallest eigenvalues of the Dirac operator are spaced AA = tt/I]V, it is 
natural to study this part of the spectrum by introducing the ruicroscopic 
variable u = VXT^, and the microscopic spectral density defined byB 

We expect that this limit converges to a universal function determined by 
the global symmetries of the QCD Dirac operator. Using partially quenched 
chiral QCD partition function we will show that Ps{u) is given by chiral Ran- 
dom Matrix Theory (chRMT) to be discussed in the next section. We will 
find the domain of validity of chRMT and will calculate the Dirac spectrum 
beyond this domain. An extensive discussion of the material of this talk with a 
complete list of references has appeared recentlyQ, and severalxxcellent recent 
reviews on Random Matrix Theory and disorder are availableu'u. 

2 Chiral Random Matrix Theory 

The chiral random matrix partition function with the global symmetries of the 
QCD partition function is defined byffl 

Z^.iM) ^ I DW jiaot (^-/ ^^y^— (3) 

where W is a nx m matrix with i> — \n~ m\ and N — n + m. As is the case in 
QCD, we assume that the equivalent of the topological charge ly does not ex- 
ceed -s/iV, so that, to a good approximation, n = N/2. Then the parameter E 
can be identified as the chiral condensate and N as the dimensionless volume of 
space time (Our units are defined such that the density of the modes N/V — 1). 
The matrix elements of W are either real (/3 = 1, chiral Gaussian Orthogo- 
nal Ensemble (chGOE)), complex {(3 = 2, chiral Gaussian Unitary Ensemble 
(chGUE)), or quaternion real {(3 = 4, chiral Gaussian Symplectic Ensemble 
(chGSE)). For QCD with three or more colors and quarks in the fundamental 
representation the matrix elements of the Dirac operator are complex and we 
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have (3 = 2. The ensembles with (3=1 and (3 = A are relevant in the case of 
two colors or adjoint fermions. The reason for choosing a Gaussian distribu- 
tion of the matrix elements is its mathematically simplicity. It can be shown 
that the correlations of the eigenvalues on the scale ofJ^hfiJmajLap^ do not 
depend on the details of the probability distributionB-ByO'EflEiy . 

3 Scales in the Dirac Spectrum 

Of course, QCD is not chiral Random Matrix Theory. This implies the exis- 
tence of a scale beyond which the correlations QCD Dirac spectra are no longer 
given by chRMT. We have studied such acaJc, ihYpJiiea of the valence quark 
mass dependence of the chiral condensate Ljlid LlEjo defined byC 
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E(m„) = (Tr — -). (4) 

niy + D 

Here, the Euclidean Dirac operator is denoted by D and the average is over the 
Euclidean QCD action which includes a fermion determinant that depends on 
the the sea-quark masses. The spectrum of D is directly related to to E(m„), 

p{\)lV=^{n^^ + e)~i:{^^-e)). (5) 

As is the case in standard chiral perturbation thcoryElH, an important scale is 
the mass where the Compton wave length of Goldstone bosons corresponding 
to the valence quark masses is equal to the linear size of the systemlHl, i.e. the 
mass scale given by 

2mcS 1 

We have identified the rangeS'B^EBS 

m^, < mc = (7) 

as the domain of validity of chRMT. In other words, below the scale rric the 
correlations of the QCD Dirac eigenvalues are completely determined by the 
global symmetries of the QCD partition function. 

A second important scale is the average position of the smallest eigenvalue 
of the Dirac operator, which is approximately given by the average spacing 
of the eigenvalues near A = 0. With spacing given by \/ p{0) and the Banks- 
Casher formula (|l]) we find 

Ami„=^. (8) 
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In the QCD Dirac spectrum we can thus distmg uishEl three different scales, 

Amin < mc < Agcn- (9) 

These scales separate four different domains. Similar domains occur in meso- 
scopic physics l3. The reason is that the transport of electrons in disordered 
samples is described by Goldstone modes. The_classical propagation of Gold- 
stone modes is given by a diffusion equation o and therefore h/nic can be 
interpreted as the time for an initially localized wave packej^to diffuse over 
the length of the sample (see for example the book by EfetovEII). The domain 
rriy <?C rric is therefore known as the ergodic domain. In mesoscopic physics the 
scale corresponding to nic is known as the Thouless energy. 



4 Partially Quenched Chiral Perturbation Theory 

The scales in QCD Dirac spectra can be studied more precisely via the generat- 
ing function of the chiral condensate. This is the partially quenched QCD par- 
tition function which, in addition to the usual quarks describing the fermion de- 
terminants, contains additional quarks and their bosonic super-partners which 
describe the queadifid. valence quark mass dependence of the chiral condensate. 
It is defined byca'cBE^ 

Nf 

where p is the Euclidean Dirac operator and Sy m [A] is the Euclidean Yang- 
Mills action. Thii. valence quark mass dependence of the chiral condensate is 
then given byE3ii3 

= Uj2 ) - T7 dj\j^,\ogZ^^imy,J). (11) 

V rriv + iAfe V 

k 

For a confining theory such as QCD, the low energy excitations of this theory 
are the Goldstone modes associated with the spontaneous breaking of chiral 
symmetry. To lowest order in the momenta, the effective action of these modes 
is completely determined by Lorentijinvariance and chiral symmetry breaking. 
In agreement with the Vafa-Wittencj theorem and maximum breaking of iihical 
symmetry E3 we expect that the chiral symmetry is broken according tooEj 

Gl{Nf + 1|1) X Gl{Nf + 1|1) ^ Gl{Nf + 1|1). (12) 

In order to obtain convergent integrals, the effective partition function is then 
based on the maximum Riemannian submanifold of the symmetric superspace 
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Gl{Nf + 1|1) which will be denoted by Gl{Nf + 1|1). It contains a kinetic 
term determined by Lorentz invariance and chiral symmetry and a mass terms 
determined by the pattern of chiral symmetry breaking. An explicit axial 
GZ(1|1) symmetry breaking is included as well 



Z{e,M) = / dUexp [ ct^x 

JueGi{Nf+i\i) J 



(13) 



In terms of the action of this partition function, the scale rric is the scale 
where the fluctuations of the mass term and the kinetic term are of equal or- 
der of magnitude. For valence quark masses <C rric the fluctuations of the 
kinetic term can be neglected and the calculation of the valence quark mass 
dependence of the chiral condensate is reduced to the calculation of zero dimen- 
sional super- integrals. Projecting onto fixed topological charge by Fourier de- 
composing the ^-dependenceiiie zero-dimensional effective partially quenched 
partition function is given bjo 

z:b{M)= f rfC/SderC/e^* ^''■(•^^+-^^"). (14) 

Ju£Gl{Nf + l\l} 

As an illustration, we mention the parameterization of the integration manifold 
for Nf = 0. Then U e Gl{l\l) is given by 

where (j) G [0, 27r], s G (— oo,cx)), and a and (3 are Grassmann variables. In 
terms of these variables the integration measure is simply given by dU = 
d(pdsdadp. The calculation of the integrals is straightforward in this case. For 
arbitrary Nf and i/, on the other hand^lhe integrals in (^) are technically 
quite involved, but we have succeeded t3E3 to obtain an analytical result for 
the valence quark mass dependence of the chiral condensate, 



^ = x{Ia{x)Ka{x) + Ia+l{x)Ka-l{x)), (16) 

where a = Nf + l^l- It coincides with the result derived from chRMT. The 
microscopic spectral density follows from the discontinuity and agrees with the 
chRMT result BQ as well. 
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Our effective partition function is amenable to cliiral Perturbation Theory 
allowing us to go beyond the Random Matrix Theory limit of QCD. Taking 
into account the kinetic term to one-loop order we have obtained an exact 
analytical expression for the QCD Dirac spectrum in the domain where chiral 
perturbation theory applies, i.e. for <C I^qcd- Among others, we have 
calculated the slope of the Dirac spectrum and shown that the spectral density 
diverges logarithmically in the quenched limit. The slope vanishes for two 
massless flavors, not only for the standard pattern of chiral symmetry breaking 
EJ, but also for other patterns of chiral symmetry breaking relevant for QCD 
with two colors and fundamental fermions or QCD with adjoint fermionsEiJ. 

The microscopic spectral density near zero virtuality was first studied in 
terms of the valence quark mass dependence of the chiral condensate. The lat- 
tice data for Yj{my) were obtained by the Columbia groupcHi for two dynamical 
flavors with sea=quark mass ma = 0.01 and TVc = 3 on a 16^ x 4 lattice. In 
Fig. 4 we plotliZi the ratio Yj[mv)l^ as a function of a; = m^VTi (the 'volume' 
V is equal to the total number of Dirac eigenvalues), and compare the results 
with the universal curves obtained from pq-QCD and chRMT (see eq. (|l6[)). 
We observe that the lattice data for different values of (3 fall on a single curve. 
Moreover, in the mesoscopic range this curve coincides with the random matrix 
prediction for Nf = v = Q. 

1.2 
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Figure 1: The valence quark mass dependence of the chiral condensate S(m) plotted as 
S(mi,)/E versus m„yE. The dots and squares represent lattice results by the Columbia 
group c3 for values of /3 as indicated in the label of the figure. 
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Of course this is no surprise. For quark masses much less than the current quark 
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mass, the fermion determinant has no bearing on the Dirac spectrum and we 
have effectively Nf = 0. At finite lattice spacing, the Kogut-Susskind action 
does not have the axial U{1) symmetry and the effects of topology become 
only visible in the continuum limit. We observe a plateau in the valence quark 
mass dependence of the chiral condensate. The increase at larger masses is 
also found in partially quenched chiral perturbation theory and is due to the 
effects of the kinetic term. The scale rric (see eq. (^)) is roughly located at the 
center of the plateau. Meanwhile, the valence quark mass dependence of the 
chiral condensate and the scalar susceptibility have also been invp s ti g ated for 
SU{2) gauge field configurations and different types of fermionsEj'Ej'EZl. In all 
cases agreement between lattice QCD spectra and the universal result (|l^) has 
been found. A detailed analysis of scalar susceptibilities in terms of quenched 
chiral perturbation theory was recently performed by Berbenni-Bitch et al. c3. 
The microscopic spectral density can be calcU|la.t.cd diri.'ptJ,A' .aiid a sin^ilar type 
of agreement with chRMT has been observedyMtHjkiyyyy ey. 
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Figure 2: The number variance T,2{n) versus n approximation for an interval starting at 
A = 0. The total number of instantons is denoted by N. 

Another interesting observable in the number variance, 'E2{n), defined as 
the variance of the number of eigenvalues in a interval containing n eigenvalues 
on average. It can be obtained in partially quenched chiral perturbation theory 
from the double discontinuity of the pseudo-scalar susceptibility cHl (A recent 
discussion of thej^lation between the pion propagator and the number variance 
can be found in P"l ) . The number variance has been calculated both by means 
lattice QCDlJo and instanton liquid E3 simulations. In Fig. 1, we show 
II3 T,2{n) versus n for eigenvalues obtained from the Dirac operator in the 
background of instanton liquid gauge field configurations. The chRMT result, 
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given by the solid curve, is reproduced up to about two level spacings. In 
units of the average level spacing, A — l/p(0) = tt/EV^, the energy Ec is given 
by ric = Ec/A = F'^L'^/tt. For an instanton liquid with instanton density 
N/V = 1 we find that Uc w 0.07%/]V. We conclude that chRMT appears to 
describe the eigenvalue correlations up to the predicted scale. 



5 Conclusions 

We have found that for mass scales where chiral perturbation theory is valid 
the fluctuations of the QCD Dirac eigenvalues can be obtained analytically 
from the partially quenched effective chiral Lagrangian. For mass scales below 
F^/SL^ the eigenvalue correlations are given by chRMT. A possible semiclas- 
sical interpretation is that quarks undergo a chaotic-diffusive motion in four 
Euclidean dimensions plus one artificial time dimension. This observation is 
consistent with a picture of confinement based on impurity scattering. 



6 Acknowledgments 

1 wish to thank the organizers of the Trento workshop and I gratefully acknowl- 
edge all my collaborators in this project. This work was partially supported 
by the US DOE grant DE-FG-88ER40388. 



References 

1. C. DeTar, Quark-gluon plasma in numerical simulations of QCD, in 
Quark gluon plasma 2, R. Hwa ed., World Scientific 1995. 

2. F.Karsch, Nucl. Phys. Proc. Suppl. 60A (1998) 169. 

3. K. Johnson, L. Lellouch and J. Polonyi, Nucl. Phys. B367 (1991) 675. 

4. T. Banks and A. Casher, Nucl. Phys. B169 (1980) 103. 

5. E.V. Shuryak and J.J.M. Verbaarschot, Nucl. Phys. A560 (1993) 306. 

6. J.J.M. Verbaarschot, Phys. Rev. Lett. 72 (1994) 2531; Phys. Lett. 
B329 (1994) 351. 



7. J.J.M. Verbaarschot, |hep-ph/9902394 



T. Guhr, A. Miiller-Groehng and H.A. Weidenmiiller, Phys. Rep. 299 
(1998) 189. 



9. G. Montambaux, |cond-mat/9602071 



10. E. Brezin, S. Hikami and A. Zee, Nucl. Phys. B464 (1996) 411. 

11. G. Akemann, P. Damgaard, U. Magnea and S. Nishigaki, Nucl. Phys. B 
487[FS] (1997) 721. 

12. G. Akemann and P.H. Damgaard, Nucl. Phys. B528 (1998) 411. 



8 



13. A.D. Jackson, M.K. Sener and J.J.M. Verbaarschot, Nucl. Phys. B479 

(1996) 707. 

14. T. Guhr and T. Wettig, Nucl. Phys. B506 (1997) 589. 

15. A.D. Jackson, M.K. Sener and J.J.M. Verbaarschot, Nucl. Phys. B506 

(1997) 612. 

16. E. Kanzieper and V. Freilikher, Phys. Rev. Lett. 78, 3806 (1997); Phys. 
Rev. E 55, 3712 (1997). 

17. J.J.M. Verbaarschot, Phys. Lett. B368 (1996) 137. 

18. J.C. Osborn and J.J.M. Verbaarschot, Phys. Rev. Lett. 81 (1998) 268. 

19. J.C. Osborn, D. Toublan and J.J.M. Verbaarschot, Nuc. Phys. B540 
(1999) 317. 

20. P.H. Damgaard, J.C. Osborn, D. Toublan and J.J.M. Verbaarschot, Nucl. 
Phys. B547 (1999) 305. ►iep-th/98H"2T^ , to be published in Nuc. Phys. 
B. 

21. D. Toublan and J.J.M. Verbaarschot, |hep-th/9904"T99| . 

22. S. Chandrasekharan, Nucl. Phys. Proc. Suppl. 42 (1995) 475; S. 
Chandrasekharan and N. Christ, Nucl. Phys. Proc. Suppl. 42 (1996) 
527; N. Christ, Lattice 1996. 

23. J. Gasser and H. Leutwyler, Phys. Lett. 188B (1987) 477. 

24. H. Leutwyler and A. Smilga, Phys. Rev. D46 (1992) 5607. 

25. B.L. Altshuler, I.Kh. Zharekeshev, S.A. Kotochigova and B.I. Shklovskii, 
Zh. Eksp. Teor. Fiz. 94 (1988) 343. 

26. A.J. McKane and M. Stone, Ann. Phys. (NY) 131 (1981) 36. 

27. K.B. Efetov, Supersymmetry in disorder and chaos, Cambridge Univer- 
sity Press, (1997); K.B. Efetov, Adv. Phys. 32 (1983) 53. 

28. A. Morel, J. Physique 48 (1987) 1111. 

29. C. Bernard and M. Golterman, Phys . Rev. D49 (1994) 486; C. Bernard 
and M. Golterman, [hep-lat/9311070| 

30. M.F.L. Golterman and K.C. Leung, |hep-lat/9711033| . 

31. C. Vafa and E. Witten, Nucl. Phys. B234 (1984) 173. 

32. M. Peskin, Nucl. Phys. B175 (1980) 197; S. Dimopoulos, Nucl. Phys. 
B168 (1980) 69; M. Vysotskii, Y. Kogan and M. Shifman, Sov. J. Nucl. 
Phys. 42 (1985) 318. 

33. J.J.M. Verbaarschot and L Zahed, Phys. Rev. Lett. 70 (1993) 3852. 

34. J.J.M. Verbaarschot, Nucl. Phys. B427 (1994) 534. 

35. M.E. Berbenni-Bitsch, M Gockeler, H. H ehl, S. Meyer, P.E. Rakow, 
A. Schafer and T. Wettig, |hep-lat/9907014. 



36. P.H. Damgaard, R.G. Edwards, U.M. Heller and R. Narayanan, hep- 



lat/9907016 



37. P. Hernandez, K. Jansen and L. Lellouch, hep-lat/9907022 



9 



38. 

39. 

40. 

41. 
42. 
43. 
44. 

45. 



46. 
47. 



48. 
49. 
50. 



M.E. Berbcnni-Bitsch, S. Meyer, A. Schafer, J.J.M. Verbaarschot and T. 
Wcttig, Phys. Rev. Lett. 80 (1998) 1146. 

P.H. Damgaard, U.M. Heller and A. Krasnitz, Phys. Lett. B445 (1999) 
366. 

M. Gockeler, H. Hehl, P.E.L. Rako w, A. Schafer and T. Wettig, Phys. 
Rev. D (in press) |hep-lat/981101g . 



R.G. Edwards, U.M. Heller and R. Narayanan, [hep-lat/9902021 
R.G. Edwards, U.M. Heller, J. Kiskis and R. Naray anan, |hcp-th/9902117 . 
F. Farchioni, I. Hip, C.B. Lang, M. Wohlgenannt, |hep-lat/9812018| . 
M.E. Berbenni-Bitsch, M. Gockeler, H. He hl, S. Meyer, P.E.L. Rakow, 
A. Schafer and T. Wettig, |hep-lat/9901013| . 

P.H. Damgaard, U.M. Heller, R. Niclasen and K. Rummukainen, hep 



lat/9907019 



F. Farchioni, I. Hip and C.B. Lang, |hep-lat/9907011 . 
M.E. Berbenni-Bitsch, M. Gockeler, T. Guhr, A.D. Jackson, J.Z. Ma, S. 
Meyer, A. Schafer, H.A. Weidenmiiller, T. Wettig and T. Wilke, Phys. 
Lett. B438 (1998) 14. 

T. Guhr, J.Z. Ma, S. Meyer and T. Wilke, |hep-lat/9806"003 . 
A. Smilga and J. Stern, Phys. Lett. B318 (1993) 531. 



A.D. Jackson and G. Carter, hep-ph/9905231 



10 



